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We study junctions between superconductors mediated by the edge states of a quantum spin Hall 
insulator. We show that such junctions exhibit a fractional Josephson effect, in which the current 
phase relation has a 4-k, rather than a 2tt periodicity. This effect is a consequence of the conservation 
of fermion parity - the number of electrons modulo 2 - in a superconducting junction, and is closely 
related to the Z2 topological structure of the quantum spin Hall insulator. Inelastic processes, which 
violate the conservation of fermion parity, lead to telegraph noise in the equilibrium supercurrent. 
We predict that the low frequency noise due these processes diverges exponentially with temperature 
r as T ^ 0. Possible experiments on HgCdTe quantum wells will be discussed. 
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Proposals for fault tolerant topological quantum com- 
putation have motivated intense current interest in find- 
ing robust physical systems that host excitations with 
non-Abclian statistics [H, [2,] . Recent experiments on the 
V = b/2 fractional quantum Hall effect have shown en- 
couraging indirect evidence for such excitations[3|,3], but 
the direct observation of non-Abelions has so far re- 
mained elusive. Recently we showed that the proxim- 
ity effect between a superconductor and a three dimen- 
sional (3D) topological insulator leads to a 2D interface 
state that supports non-Abelian Majorana fermions 5]. 
A first step towards implementing this proposal would 
be to demonstrate experimentally the topological order 
responsible for Majorana fermions. 

In this paper we study Josephson junctions mediated 
by a 2D topological insulator, also known as a quan- 
tum spin Hall insulator (QSHI) [6, 7, 8, 9, 10]. We pre- 
dict that such junctions exhibit a fractional Josephson 
effect^ whose origin is related to the presence of Majo- 
rana fermions. The signature of the fractional Josephson 
effect is that the current phase relation has a An rather 
than a 27r periodicity. This behavior was first predicted 
by Kitaev using an idealized model of a ID spinless p 
wave superconductor [l3|. Kwon et al.[l3] proposed that 
a related effect can occur at junctions between uncon- 
ventional 3D superconductors. They argued that it leads 
to an AC Josephson effect with half the usual Josephson 
frequency and that in a weak tunneling limit the Joseph- 
son current is carried by electrons rather than Cooper 
pairs. The 47r periodicity can occur because the junction 
has two states with different Josephson currents that are 
interchanged when the phase is advanced by 27r. At finite 
temperature inelastic processes can cause transitions be- 
tween the states, leading to telegraph noise in the Joseph- 
son current. We will show that in our setup (unlike Ref. 



12f ) these transitions are forbidden by the local conserva- 



tion of fermion parity, which counts the number of elec- 
trons modulo 2. This leads to an exponential suppression 
of the transition rate at low temperature. This can be 
experimentally probed by measuring the low frequency 



current noise S{ijJ 0), which we predict diverges expo- 
nentially at low temperature. 

The QSHI is a time reversal invariant insulating 
state with a bulk energy gap generated by spin or- 
bit interactions 0, It has recently been observed 
in HgCdTe quantum wells [l^. The QSHI is distin- 
guished from an ordinary insulator by a Z2 topological 
invariant , which necessitates the existence of gapless 
edge states. The states at the edge of the QSHI form a 
unique ID system that is essentially half of an ordinary 
spin degenerate ID electron gas. In the simplest case it 
consists of a single band of right moving electrons paired 
via Kramers theorem with a left moving band with the 
opposite spin. These states are robust in the presence of 
disorder because time reversal symmetry prevents elastic 
backscattering. In the absence of inelastic scattering the 
transmission in the edge states is perfect. 

Suppose the edge states are in intimate contact with 
an s wave superconductor. Due to the proximity effect, 
the tunneling of Cooper pairs will induce a pairing poten- 
tial A = Aqc*"^ in the edge states, which depends on the 
phase (j) of the superconductor and the nature of the con- 
tact. Using the notation of Ref. 5' write H = '^^H'^ /2, 
where ^ = [(-01, i/ii), (?/;|, — i/jI)] is express in terms of 
field operators V'T(i) describing the right(lcft) movers and 



H = -ivT^a^dx " ti-Tz + Ao(cos(/)T2; 



- sin (j)Ty ) . 



(1) 



aj are Pauli matrices acting in the space of right and left 
movers '(/'T.i ^-nd Tj are Pauli matrices which mix the ^ 
and ■0^ blocks of 4". u is the velocity of the edge states, 
fi is the chemical potential and we set h = 1. The eigen- 
states of ([1]) come in pairs at ztE. Due to the redundancy 
in 5", these states are not independent, and the Bogoli- 
ubov quasiparticle operators satisfy F-b = F^. 

Eq. [1] is similar to Kitaev's model of superconducting 
spinless electrons in 1D[iI[. In Kitaev's model there are 
zero energy Majorana bound states associated with the 
ends of the sample. In our system, the edge - which is 
the boundary of the 2D QSHI - can not have an end. 
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FIG. 1: A S-QSHI-S junction in an RF SQUID geometry 
where the QSHI forms a Corbino disk. 



By breaking time reversal symmetry, however, a Zeeman 
field can introduce a mass term into H of the form 



(2) 



When M > Vz opens an insulating gap in the edge 
state spectrum. Vz could arise either from an applied 
magnetic field (as in Ref. [13]) or due to proximity to a 
magnetic material. Zero energy Majorana bound states 
will exist at the interface between regions with gaps dom- 
inated by A and mQ. In the presence of both A and 
M the gap is the smaller of |Ao ± M\. When Aq = |Af | 
a single band is gapless, and for Aq ~ |M| the low en- 
ergy sector of ([T]) has the form of a Su-SchriefFer-Heeger 
model[l3|. which has a well known zero energy bound 
state where Aq — \M\ changes sign. The Bogoliubov 
quasiparticle operator associated with this state is a Ma- 
jorana fermion, which satisfies 70 = 70 . 

Consider a superconductor-QSHI-superconductor (S- 
QSHI-S) junction in which the edge states of a QSHI 
connect two superconductors separated by a distance L. 
Fig. 1 shows an RF SQUID geometry, in which the phase 
difference across the junction (j) = {2e/K)^ is controlled 
by the magnetic flux We also assume that the QSHI 
forms a Corbino disk which circles the flux. As a prac- 
tical matter, this geometry is not essential, but we will 
see that it has considerable conceptual value. We will 
also consider the effect of a Zeeman term in the gap be- 
tween the superconductors, which will make the connec- 
tion with Majorana bound states transparent. We em- 
phasize, however, that there will be a non trivial effect 
even when this term is absent. To determine the char- 
acteristics of the junction we solve the Bogoliubov de 
Gennes (BdG) equation with 



A(a;) = Ao [6'(-a;-L/2) + e"^6'(x-L/2)] 
M{x) = MQe{x + Lli)6{-x + L/2). 



(3) 



By matching the solutions it is straightforward to de- 
termine the spectrum of Andreev bound states in the 




FIG. 2: Spectrum of Andreev bound states in the junction as 
a function of phase difference <j) for parameters indicated in 
each panel. L is in units of «/Ao and Mo and ^ are in units 
of Aq. (a) and (c) are independent of 



junction. The calculation is similar to Ref. [l^i as well 
as the theory of superconducting quantum point con- 
tacts(SQPCs) [3, [IB]. However, we shall see that there 
is a fundamental difference with both of those theories. 

Fig. 2a shows the spectrum as a function of for A/q — 
0. For L < f/Ao there is a single pair of bound states 
E = ±eo ((/)). For L <C w/Aq our model reduces to the 5 
function model solved in Ref. [12, where the normal state 
transmission probability is D — l/[l-|-(Mosinh(KL)/K)^], 



with K = \J Mq — /i^. In that case 



eo(</') = Vi:'Aocos(0/2). 



(4) 



Fig. 2b shows a case where Mq ^ Aq, so the normal 
state transmission D < 1. When D <C 1 there are two 
weakly coupled Majorana end states at a; = ±L/2. When 
L > v/ Aq there will be additional Andreev bound states 
in the junction with a level spacing of order v/L. Fig. 2c 
shows the case where vL/Aq — 3 with Mq — 0, in which 
time reversal symmetry requires Kramers degeneracies 
when (j> — or IT. Fig. 2d shows the effect of finite Mq 
and /i, which lifts most of the degeneracies. However, the 
crossing at -E = remains, and is of special significance. 

To understand the level crossing consider a low energy 
theory for _B <C Ag. The two eigenvectors of (HJSl) 
with energy ±eo((^) define Bogoliubov operators Foi ~ 
^"^"^oi- Due to particle-hole symmetry, Fo+ = rj_ = Fq. 
The low energy Hamiltonian is thus 



H = eo((/))(r5ro - 1/2) = 2ieo((/.)7i72, 



(5) 



where we introduced Majorana operators 71 = {T+r^)/2, 
72 = -«(F - r'f)/2. In the weak tunneling {D < 1) 
limit 7i_2 describe the Majorana end states at a; = ±L/2, 
which are weakly coupled by tunneling of electrons. Eq. 
^describes two states distinguished by A^o = tJFo = 0, 1. 
Mixing these states requires an interaction that changes 
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Nq. Due to the pairing term in ([T]), the total charge is not 
conserved. However, the fermion parity, defined as the 
number of electrons modulo 2, is conserved in ([T][3]). This 
forbids the coupling between the two states and protects 
the crossing at eo(0) = 0. 

There is a problem, however, with the fermion parity. 
The junction Hamiltonian ((T][3]) is invariant under a 2tt 
phase change, but when — > ^ + 27r, the system passes 
through a single level crossing and can only return to 
the initial state by a process which changes Nq by 1. 
The fermion parity thus apparently changes when (p — > 
+ 2n. This has to do with the unbounded spectrum as 
E ^ —oo and reflects a, fermion parity anomaly similar to 
the SU{2) anomaly in 4D field theory |16|]. This anomaly 
is related to non-Abelian statistics. When (j) advances by 
27r 7i 7i and 72 — > —72. In the tunneling limit this 
can be interpreted as Ivanov's rule[l3| for the effect of 
braiding a vortex between the Majorana bound states. 

The physical origin of the fermion parity anomaly lies 
in the topological structure of the QSHI. Consider first 
the Corbino disk in Fig. 1 without the superconductor. 
In Ref . Ts' we showed that the Z2 invariant characterizing 
the QSHI describes the change in the Z2 "time reversal 
polarization" (TRP) when flux h/2e is threaded through 
the hole. A nonzero TRP specifles a many body Kramers 
degeneracy localized at either edge of the disk. Since an 
odd number of fermions has a Kramers degeneracy, the 
TRP is precisely the fermion parity. With the super- 
conductor present, start in the groundstate at $ = 0. 
When flux <i> = ft,/2e is threaded through the hole, (f) ad- 
vances by 27r and a unit of fermion parity is transferred 
from the inner edge of the disk to the junction on the 
outer edge. The anomaly in (1-3) is similar to the chi- 
ral anomaly for edge states in the quantum Hall effect, 
where bulk currents violate charge conservation at the 
edge. Note however that though (1-3) is invariant under 
(j) ^ + 27r, the global Hamiltonian, which includes the 
bulk QSHI, is physically distinct when $ = and h/2e. 

The local conservation of fermion parity has important 
consequences for the current and noise in a S-QSHI-S 
junction. This is most striking near the degeneracy point 
for eo ^ Ao and T <C Aq. For the remainder of the paper 
we will focus on that regime. We will also consider the 
limit L <^ v/Aq, where there is a single Andreev bound 
state and (|4]) applies, though the results can straightfor- 
wardly be generalized to the case with multiple Andreev 
levels, provided T <^v/L. In this case, A^o distinguishes 
two states, with Josephson currents I± = ±Iq with 



/o(0) = i/DeAosin0/2. 



(6) 



In the absence of transitions that violate local fermion 
parity conservation there can be no transitions between 
/+ and /-, signaling a fractional Josephson effect. 

Elastic scattering processes can be incorporated into 
the BdG Hamiltonian from the start, and will not lead to 



violations of the fermion parity. However, at finite tem- 
perature, inelastic processes j2li [2^ can lead to a transi- 
tion between and provided an available fermion 
is present to switch the fermion parity. This could be 
either due to a thermally excited quasiparticle or due to 
hopping from a bulk localized state. These processes, 
however, will be exponentially suppressed at low tem- 
perature. On a time scale longer than the switching time 
the current will thermalize, with an average value [il 11^ 



(/(0)) =/o(0)tanheo(0)/2T. 



(7) 



On shorter times, the current will exhibit telegraph noise, 
as it switches between I±. 

In order to model the inelastic processes responsible 
for the telegraph noise we consider the interaction of the 
Andreev level Fq with a bath of fermions c„ (e.g. quasi- 
particles) and bosons bm (e.g. phonons). We thus write 



+ VnmCnbljTQ + h.c] 



(8) 



Here En, a;„ > 0, and we have ignored terms which create 
(or anihilate) both fermions and bosons. The transition 
rates r^^ (cq, T) between the states Nq and Nq ± 1 follow 
from Fermi's golden rule. For eo,T <C Aq we find 



.-1 



,± -eT'«/2^ (i«i(T)e^<'/2?^ + u.2(r)e-'«/2^) , (9) 
where 

w^MT) - 2nY,e-^-/^\V^i\'SiEn - c^™). (10) 

If either the Zeeman term vanishes {Mq = 0) or the sys- 
tem is symmetric under x —^ —x, then wi(T) = W2{T) = 
w(r). We will assume this below, though the results are 
only slightly modified otherwise. w{T) depends on the 
dominant source of fermions, which we take to be either 
thermally activated quasiparticles or Mott variable range 
hopping from bulk localized states. 



-Aa/T 



w{T) cx ';_(T„/T)V3 



quasiparticles, 
hopping. 



(11) 



To depends on the density of states and localization 
length, and we assume the hopping is 2D. 

The transition rate is exponentially suppressed for 
T ^ 0. At sufficiently low temperature the result- 
ing telegraph noise could be observed in the time do- 
main. At higher temperature there is a signature in 
the noise spectrum S{lu). We determine S(uj) semiclas- 
sically by solving a kinetic equation for the probability 
p{t) that A^o = ijlll liil. This has the form dp/dt = 



-{p — p)/t, where t ^ = 



Aw cosh^ eo/2T. 



P 



(1 -I- expeo/T) ^ follows from the detailed balance 
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condition t^/t- = e'^"/'^ . Temporal correlations in I{t) 
decay exponentially on a time scale w"^, and the noise 



spectrum S{uj) = 2 e*"*(/(t)/(0)) is given by[ig|,l2l| 



4/2 



cosh eo(0)/2rl + w2r2 
In the zero frequency limit we have 



S{uj 0) 



-'o 



w(T)cosh^eo(0)/2T 



(12) 



(13) 



For D = 1, these results are similar to the theory of 
a SQPC[il, [23, mi, 111. However, the current in Eq. 



([6]) is half the value of a perfect single channel SQPC. A 
SQPC is like two copies of a S-QSHI-S junction. This 
leads to a fundamental difference, because in the SQPC 
there is no conservation law to prevent scattering between 
the ±/o states, which can occur via low energy processes 
that transfer an electron between the two pairs. Elastic 
backscattering in the SQPC leads to an avoided crossing 
of the states near = 0, and the Andreev states carry no 
current at (p = tt. Even if the transmission of the SQPC 
is perfect, inelastic processes will couple the states. For 
instance, near the degeneracy point eo((/') = spin flip 
scattering via the nuclear hyperfine interaction will lead 
to a finite lifetime for transitions, even at low tempera- 
ture. It is also of interest to compare with the theory of 
Ref. [1^ In that work, multi channel junctions were con- 
sidered. Independence of the different channels requires 
translational symmetry, so impurity scattering will lead 
to the violation of the conservation fermion parity within 
a given channel. Thus the low temperature behavior pre- 
dicted by ([TT|) and is unique to the S-QSHI-S junc- 
tion, and is a signature of the fermion parity anomaly. 

We now briefly consider junctions at finite voltage bias. 
The analysis is similar to Ref. [3. There are two cases 
to consider, depending on Mq. For Mq = 0, the perfect 
transmission of the edge states causes the Andreev levels 
to merge with the continuum levels. This leads to a finite 
DC current, which for eV <^ Aq can be understood semi- 
classicalUy in terms of multiple Andreev reflections. For 
w{T) < el/ < Ao, the current is I{V) = {2/Tr)IcSgnV, 
where Ic = VDeAo/2. For Mq > 0, there is an en- 
ergy gap S separating the Andreev levels from the con- 
tinuum, as in Fig. 2(b,d). For w{T) eV S there 
will be a fractional AC Josephson current with frequency 
ey/?i[l^. For eV ~ S Landau-Zener tunneling processes 
through the gap S will lead to a damping of the AC 
Josephson current as well as a finite DC current. 

We close by discussing the feasibility of experiments 
using the QSHI recently achieved in HgCdTe quantum 
wells which has a bulk gap of order 20 meV 23|. 

The desired geometry would be similar to Ref. [2J], 
where a 2D InAs quantum well was contacted with Nb. 
The proximity induced gap will depend on the con- 
tact. If optimized Aq could be of order the bulk gap 



of the superconductor. To determine the required junc- 
tion size we use v = 3.6 eVA 23| and Aq = .ImeV. 
Then, L < v/Aq ~ 3/im sets the scale for having a sin- 
gle Andreev level. The simplest experiment would be 
to study a single current biased junction, which is pre- 
dicted to have a critical current Ic ~ eAo/2 ^ lOnA, 
which is half the value of a perfect single channel SQPC. 
Measuring the equilibrium noise at (/> tt requires an 
inductive measurement on a ring[l3|. The physics at 
Mq =/= requires a magnetic field in the junction region. 
For an appropriately aligned field, a field induces a gap 
B X (3.1meV/T)[2^ in the edge states, so a field of order 
.03r could suppress the normal state transmition D as 
well as the magnitude of the Josephson current. 
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